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EXECUTIVE  SUMMARY 


Press  reports  of  "East  Coast  Mystery  Booms"  have  led  to  a number  of 
studies  of  the  propagation  of  shock  waves  into  the  thermosphere.  The 
energy,  E , in  an  N-wave  is  proportional  to  the  square  of  the  relative 
overpressure,  s , 

v , 2 

E - Lps 
s » Ap/p 

where  L is  the  length  of  the  N-wave,  p the  ambient  pressure  and  £p 
the  shock  overpressure.  The  fractional  energy  loss  due  to  the  entropy 
jump  across  the  shock  traveling  in  the  x direction  is  proportional  to 
the  relative  overpressure  at  every  point  along  the  travel  path 

1 dE  s 

E dx  2,'L 

where  y is  the  ratio  of  the  heat  capacities.  Solving  for  the  relative 
overpressure  at  high  altitudes,  we  obtain  a limiting  strength  of  the  shock 
of 

3L 

H 


where  H is  the  scale  height.  Even  though  the  N-wave  lengthens  as  the 
shock  propagates  upwards,  the  increase  in  scale  height  insures  that  shock 
remains  weak  at  all  heights.  For  a typical  Concorde  Mach  2 boom  having 


ill 


an  initial  strength  of  100  Pa,  the  turnover  height  of  the  boom  is  160  km. 

In  reaching  this  altitude,  902  of  the  energy  is  lost  by  100  km  and  97%  of 

the  energy  is  lost  by  160  km.  As  a result  of  the  shock  remaining  weak, 

_2 

the  change  of  ambient  temperature  is  only  about  10  0K  and  the  velocity 
imparted  to  the  thermosphere  is  about  3 cm/s.  The  boom  is  also  lengthened 
to  about  4 km,  so  that  the  peak  frequency  at  ground  levels  is  about  0.1  Hz. 
SST  generated  booms  thus 

1)  do  not  heat  the  thermosphere; 

2)  do  not  create  thermospheric  wind; 

3)  do  not  cause  chemical  changes  in  the  thermosphere. 

General  considerations  of  the  propagation  of  infrasound  suggest  a 
number  of  research  priorities  including: 

1)  construction  and  operation  of  directional  infrasound 
arrays  to  monitor  natural  and  artificial  infrasound 
sources; 

2)  measurement  of  the  source  spectrum  of  a wide  variety 
of  artificial  sources; 

3)  coupled  lnfrasonic  and  ionospheric  observation  to 
study  the  interaction  of  near  surface  generated 
disturbances  with  the  ionosphere; 

4)  construction  and  operation  of  a directional  array  of 
lnfrasonic  sources  to  study  propagation  conditions 
in  the  thermosphere. 
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OVERVIEW 


Public  attention  was  drawn  to  unexplained  atmospheric  sonic  phenomena 

by  a series  of  press  reports  of  "Mysterious  East  Coast  Booms"  from  early 

1 2 

December  1977  through  March  1978  * . Numerous  citizen  reports  of  intense 

sounds  in  New  Jersey,  in  the  vicinity  of  Charleston,  South  Carolina  and 

on  the  southwest  tip  of  Nova  Scotia  prompted  a number  of  independent 
345 

investigations  * * . In  addition  to  the  citizen  reports,  certain  of  the 
events  were  recorded  by  acoustic  arrays  of  the  Lamont -Doherty  Observatory^ 
and  by  seismic  arrays  of  Lamont-Doherty  Observatory,  Western  Observatory 
and  Baptist  College  in  Charleston  where  acoustic  events  are  distinguished 
by  vertical  motion  with  no  horizontal  displacement.  Citizens  also  noted 
lights  in  the  sky,  sometimes  associated  with  the  noises  at  times  when 
there  was  no  thunderstorm  activity.  Sound  events  have  also  been  reported 
in  other  areas  of  the  east  coast  such  as  the  Tidewater  region  of  Virginia, 
but  almost  all  reports  have  come  from  Nova  Scotia,  New  Jersey  and  South 
Carolina.  Press  reports  of  events  in  New  Jersey  and  South  Carolina  ceased 
in  late  March,  but  citizen  reports  of  booms  in  Nova  Scotia  continued 
through  August  1978. 

Public  concern  led  Senator  Harrison  Williams  (N.J.)  to  request  the 
White  House  to  conduct  an  investigation.  Frank  Press,  Director  of  the 
Office  of  Science  and  Technology  Policy,  referred  the  matter  to  the 
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Department  of  Defense  which  In  turn  designated  the  Naval  Research  Labor- 

3 

atory  (NRL)  as  the  lead  agency.  The  resulting  NRL  study  concluded  that 
a majority  of  the  New  Jersey  and  South  Carolina  booms  could  be  explained 
by  the  supersonic  flight  of  military  aircraft  in  the  "warning  areas"  off 
the  east  coast.  The  NRL  study  noted  that  unusual  temperature  and  wind 
conditions  could  probably  account  for  the  amplification  and  propagation 
of  the  sonic  carpet  boom  over  the  long  distances  required  if  the  military 
aircraft  were  responsible.  The  NRL  study  also  briefly  noted  that  the 
booms  in  Nova  Scotia  could  be  attributed  to  Concorde  flights. 

In  its  study,  NRL  examined  many  alternative  explanations  for  the  booms, 
including:  nuclear  explosions,  military  research  and  development  activities, 
military  or  civilian  use  of  high  explosives,  ship  disasters,  USSR  ship  oper- 
ations, geophysical  exploration,  antipodal  events,  missile  launches  and  re- 
entries, low  altitude  satellites,  explosive  freezing  of  supercooled  water 
vapor,  meteorites,  winter  lightning,  direct  seismic  generation  and  the 
explosion  of  biogenic  or  abiogenic  methane.  NRL  concluded  that  all  man- 
made sources  could  be  excluded  with  the  exception  of  military  aircraft  off 
Nova  Scotia.  NRL  considered  that  natural  sources  were  highly  improbable; 
a conclusion  that  has  been  disputed  by  Chalton  and  MacDonald  in  a compre- 
hensive  analysis  of  about  600  separate  events  . 

Noting  that  the  reports  of  the  east  coast  booms  approximately  coin- 
cided with  the  Initiation  of  Concorde  service  into  New  York's  J.F.  Kennedy 
Airport  in  late  November  1977,  Jeremy  Stone  of  the  Federation  of  American 
Sclentista  suggested  that  Concorde  was  responsible  for  the  New  Jersey  and 
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Charleston  events  as  well  as  for  the  Nova  Scotia  booms.  If  Concorde  were 
the  source  of  the  South  Carolina  booms,  then  the  shock  would  have  to  propa- 
gate 1200  km,  the  distance  from  the  last  supersonic  segment  of  the  flight. 
Liszka^  has  observed  at  Kiruna,  Sweden,  Concorde  flights  at  distances  of 
3000  to  4500  km  with  a rms  amplitude  of  about  0.1  Pa  in  the  frequency 
band  around  2 Hz.  Liszka  interprets  the  observed  signal  as  a focused 
boom  in  which  the  signals  from  a finite  segment  of  the  flight  path  arrive 
simultaneously  at  a particular  point  on  the  earth's  surface. 


Stimulated  by  Stone's  speculation,  Garwin  carried  out  a preliminary 
analysis  of  the  propagation  of  a linear,  lossless  wave  into  the  tenuous 
and  hot  thermosphere.  A Concorde  sonic  boom  has  a leading  shock  with  an 
overpressure  Ap  of  about  100  Pa  near  ground  level;  that  is,  the  strength 

_3 

(relative  overpressure)  which  we  designate  by  s * Ap/p  , is  about  10 
With  Garwin 's  approximations,  such  a shock  becomes  very  strong  at  180  km 
altitude,  since: 


S180 

s 

g 


7 * 103 


where  p is  the  ambient  pressure  and  where  the  energy  E in  the  shock 

2 

is  proportional  to  (Ap/p)  . The  corresponding  temperature  rise  inside 

4 

the  shock  would  be  10  #K  with  a resulting  kinetic  temperature  correspon- 
ding to  a particle  velocity  well  above  the  escape  velocity.  Garwin  tenta- 
tively concluded  that  the  Concorde  generated  shocks  could  significantly 
perturb  the  thermosphere  by  depositing  energy  and  momentum  leading  to 
thermospheric  escape  and  high  thermospheric  winds.  Further,  the  shock 


would  propagate  through  the  thermosphere  at  an  average  velocity  slightly 
less  than  the  velocity  of  the  aircraft  at  the  point  from  which  the  shock 
was  launched.  Thus  the  acoustic  disturbance  could  travel  through  the 
atmosphere  at  mean  velocities  well  above  the  velocity  of  sound  at  sea 
level.  The  booms  refracted  by  the  thermosphere,  called  hyperbooms  by 
Garwin,  could  travel  horizontally  for  distances  of  hundreds  of  thousands 
of  kilometers  and  through  focusing  might  produce  the  observed  east  coast 
booms . 

9 

Gardner  and  Rogers  have  carried  out  a detailed  analysis  of  the  propa- 
gation of  a weak  shock  in  the  atmosphere  taking  explicit  account  of  the 
non-linear  character  of  even  a weak  shock,  the  density  and  thermal  struc- 
ture of  the  atmosphere  and  attenuation  due  to  linear  dissipation.  Gardner 
and  Rogers  find  that  the  non-linearities  of  a weak  shock  lead  to  a substan- 
tial energy  loss  as  the  shock  propagates  upward  in  the  thermally  stratified 
atmosphere.  They  conclude  that: 

1.  Ap/p  never  exceeds  0.2  in  the  thermosphere  so  that  the 
shock  remains  weak  throughout  its  path.  A weak  shock 
could  neither  significantly  perturb  the  temperature  or 
wind  field  of  the  thermosphere. 

2.  A Concorde  boom  generated  at  17.5  km  height  with  the 
airplane  traveling  at  Mach  2 will  reach  the  ground  at 
a horizontal  distance  of  320  km  from  the  airplane. 

3.  The  sound  received  at  this  distance  will  have  a small 
overpressure  of  less  than  0.5  Pa  and  the  energy  will 
peak  at  about  0.2  Hz. 


A 


» 


The  analysis  present  d in  this  report  carried  out  independently  of 
the  NRL  work  uses  a much  simpler  model  of  sonic  boom  propagation  to  ar- 
rive at  conclusions  similar  to  those  of  Gardner  and  Rogers.  In  this 
analysis  the  fractional  energy  loss  due  to  the  entropy  jump  for  a shock 
propagating  in  the  x direction  is 


I 


I 


t 


1 dE  ~ (Ap/p) 

E dx  2 y L 

where  L is  the  sonic  boom  length  and  y the  ratio  of  the  heat  capaci- 
ties. The  fractional  energy  loss  is  thus  proportional  to  the  overpressure 
at  every  point  along  the  travel  path.  The  relative  overpressure  Ap/p 

stays  well  below  the  limit  of  2yL/H  at  all  altitudes,  where  H is 

s s 

the  scale  height  of  the  atmosphere.  This  limit  stays  well  below  0.5  at 
all  altitudes  (even  though  the  sonic  boom  lengthens  as  it  travels  to  high 
altitudes)  because  the  scale  height  of  the  atmosphere  also  increases  at 
high  altitudes.  Thus,  the  large  scale  height  in  the  thermosphere  insures 
that  the  shock  will  remain  small. 


The  results  obtained  here  and  by  Gardner  and  Rogers  are  consistent 
® 10 

with  Grover's  observations  of  Concorde  generated  shocks.  Using  an  infra- 
sonic  array,  Grover  measured  at  a horizontal  distance  of  300  km  from  the 
Concorde  flight  path  a pressure  pulse  of  0.4  to  0.8  Pa  in  the  frequency 
range  of  0.3  to  0.5  Hz. 
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Liszka's  observations  of  pressure  amplitude  of  0.1  Pa  at  a distance 

of  4500  km  and  at  a frequency  of  about  2 Hz  can  also  be  interpreted  in 

terms  of  the  calculations  presented  here.  After  the  first  bounce  off  the 

thermosphere,  the  ground  pressure  (Ap/p)  at  a horizontal  distance  of 

300  km  is  about  10  ^ . The  ground-reflected  shock  strength  is  smaller  by 

a factor  of  100  compared  to  the  initial  shock  having  a pressure  amplitude 

of  100  Pa.  With  the  small  over-pressure,  non-linear  weakening  will  be 

small  and  the  infrasound  wave  will  propagate  a distance  R with  only 

-U 

geometrical  reduction  in  strength  proportional  to  R provided  the 
attenuation  is  negligible.  For  infrasonic  waves  attenuation  at  heights 
below  100  km  is  less  than  about  .001  db/kn  at  1 Hz^.  For  distances  of 
4500  km,  the  geometrical  spreading  reduces  the  pressure  amplitude  by  a 
factor  of  four  and  focusing  can  be  called  on  to  overcome  dissipation 
losses  above  100  km.  Thus  a measured  rms  amplitude  of  0.1  Pa  at  4500  km 
distance  is  consistent  with  an  amplitude  of  0.4  to  0.8  Pa  at  300  km  even 
though  the  longer  propagation  involves  some  14  bounces  off  the  thermosphere. 
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II 


PROPAGATION  OF  SONIC  BOOMS  IN  A STRATIFIED  ATMOSPHERE 


» 
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We  show  here  that  upward-going  sonic  booms  will  not  significantly 
affect  thermospheric  temperatures  or  winds,  and  that,  although  these 
upward-going  booms  may  be  refracted  downward  by  the  thermosphere,  they 
lose  so  much  energy  during  their  upward  travel  that  they  cannot  result  in 
large  sonic  booms  on  the  ground. 

For  simplicity  of  illustration  we  consider  vertically  traveling 
planar  sonic  booms.  We  thus  avoid  the  necessity  for  considering  either 
geometrical  spreading  or  the  change  in  direction  due  to  refraction  that 
a non-vertical  wavefront  would  experience.  We  can  draw  all  significant 
conclusions  from  this  simple  case.  Our  basic  problem  is  to  solve  for  the 
sonic  boom  strength  as  a function  of  altitude.  We  know  that  rather  close 
to  the  aircraft  altitude,  the  strength  (relative  overpressure)  is  about 

_3 

10  . If  the  boom  reaches  a strength  of  unity,  a number  of  important 

effects  would  occur;  but  if  our  vertically  traveling  sonic  boom  always 
has  a strength  well  below  unity  as  it  travels  into  the  rarefied  upper 
atmosphere,  then  no  sonic  boom  at  any  angle  will  ever  become  strong 
(ignoring  cusps  and  caustics),  and  we  can  easily  calculate  sonic  boom 
effects. 
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l.  Sonic  Boom  Characteristics 

A supersonic  aircraft  creates  a complicated  shock-wave  structure  that 
propagates  away  from  the  aircraft's  flight  path.  At  moderate  distances 
from  the  aircraft,  the  structure  is  reasonably  well  approximated  by  a 
symmetric  N-wave  as  shown  in  Fig.  1. 


Pressure 


Figure  1 PROPAGATION  OF  AN  IDEAL  N-WAVE 


* The  true  sonic  boom  will  differ  from  this  idealized  form,  but  the 

differences  will  not  affect  our  results  (see  Appendix  A). 


i 


Once  we  have  established  the  pressure  waveform,  we  may  calculate 
various  useful  quantities  from  elementary  Bhock-wave  theory  (see 
Appendix  A) : 


n 
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h « altitude  (distance  along  direction  of  sonic  boom  travel) 
u * speed  of  first  shock  discontinuity 

c ■ speed  of  sound  * (yR^T) 

= ambient  temperature  before  sonic  boom  passage 
Pa  « ambient  pressure  before  and  after  sonic  boom  passage 
Ap  = overpressure  just  behind  the  first  shock  discontinuity 
s - relative  overpressure  Ap/p  (called  the  sonic  boom  strength) 

3 

y = ratio  of  specific  heats  for  the  atmosphere  as  an  ideal  gas 

AT  = increase  in  temperature  just  behind  the  first  shock 

discontinuity 

AT  = increase  in  ambient  temperature  after  the  passage  of 
the  entire  sonic  boom 

E « total  energy  (per  unit  area)  in  the  sonic  boom  wave  front. 
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Within  the  boom  the  temperature  rise  is  first  order  in  the  strength 


AT  __  y-1 
T * Y 


(1) 


while  the  heating  of  the  air  behind  the  boom  is  third  order: 


ATa  y 2-l  s3 
Ta  12?- 


(2) 


The  energy  in  the  sonic  boom  is  given  by 


E L p s2 

6y^  a 


(3) 
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As  the  sonic  boom  travels,  energy  is  lost  in  heating  the  atmosphere.  The 


i 

I 


energy  loss  relation  is 


dE  ~ s 
E dh  -2yL 


(4) 


or  twice  as  large  if  the  ideal  N wave  shape  is  maintained  (see  App.  A). 


The  air  behind  the  spnic  boom  will  not  only  be  heated  (by  AT  ) 
but  will  also  be  given  momentum.  This  momentum  transfer  will  result 
in  atmospheric  winds,  with  speeds  given  by 


AT 

_ c a 

v “ y(y-i!)  t 

a 


(5) 


The  first  shock  travels  slightly  faster  than  an  acoustic  wave: 


— / Y+l  \ 

u"c(l  + 17s) 


(6) 


resulting  in  a lengthening  of  the  boom  while  it  travels 


dL  _ u-c  — Y+l 
dh  “ c "" 


(7) 


The  above  simple  relations,  plus  a model  for  the  atmosphere,  are 
all  that  is  needed  to  establish  the  results  of  this  section. 
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B.  No  Lengthening;  Isothermal  Atmosphere 

A simple  example  will  Illustrate  the  physics  of  the  situation. 

Let  us  ignore  the  changes  in  L with  distance,  taking  L as  a constant 
Lq  . Furthermore,  let  us  consider  an  isothermal  atmosphere.  The  sound 
speed  cq  is  a constant  in  such  an  atmosphere.  We  need  the  pressure  as 
a function  of  altitude  in  order  to  make  use  of  (3).  We  obtain  this 
through  the  sound-speed  relation 


t 


t 


t 


c2  - YP/P 


and  the  hydrostatic  relation 


3hP  * ~Pg 

resulting  in  the  solution 


p 

a 


-h/H 

s 

V 


(8) 


(9) 


(10) 
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That  is,  the  well-known  exponential  atmosphere.  Let  Pq  , be  the 
ambient  pressure  and  sonic-boom  strength  at  the  aircraft  altitude. 

Then  (3),  (4),  and  (10)  yield  a differential  equation  for  the  energy  E 
as  a function  of  altitude: 


The  solution  to  this  equation  can  be  easily  found,  and  yields  immediately 
the  energy  and  strength  as  a function  of  altitude: 


H 

/h/2Hs  M 

E * E 

o 

1 + TT~  s 
2yLo  o 

k i 

h/2H 

s * 

s e 
o 

H ! 

1 + s 

2yL  o 

o 

" h/2H 

s i 

e -1 

The  solution  for  the  strength  shows  a strong  saturation  phenomenon, 
since  at  high  altitude,  regardless  of  the  initial  strength  at  aircraft 
altitude,  we  have 

s_!^o  = 

h-+«  H ' max 
s 

Typical  values  for  L (~300  meters)  and  H (~10  km)  will  yield  a 

o s 

saturation  value  s ~0.1. 

max 


Ue  could  at  this  point,  easily  calculate  all  the  atmospheric  effects 
of  interest  within  this  simple  model.  We  would  find  that  a strength  of 
<0.1  never  leads  to  significant  heating  or  winds  in  the  upper  atmosphere 

_3 

Eqs . (1),  (2)  and  (5),  and  that  a boom  with  sQ~10  loses  more  than  95% 
of  its  energy  by  the  time  it  has  reached  a height  of  140  km  Eq.  (12). 
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(14) 
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C.  Lengthening 

However,  the  lengthening  of  the  boom  Eq.  (7)  is  not  negligible.  If 
we  allow  the  boom  to  lengthen  in  our  calculations,  then  as  the  altitude 
increases  the  saturation  value  of  the  strength  also  increases.  As  a 
result  the  boom  strength  may  continually  increase;  with  typical  initial 
parameters  we  find  that  the  boom  would  become  a strong  shock  (s  — 1)  at 
170  km  altitude. 


t 


t 


t 


t 
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The  increase  to  a strong  shock  occurs  because  of  the  rapid  decrease 
in  ambient  pressure  in  the  assumed  isothermal  atmosphere.  We  find  that 
if  a more  accurate  model  of  the  atmosphere  that  includes  the  thermosphere 
is  used,  the  boom  does  not  become  strong. 

D.  Atmospheric  Model 

The  temperature  in  the  atmosphere  (Fig.  2a)  begins  to  increase 
rather  suddenly  at  an  altitude  of  about  100  km.  We  model  this  increase 
by  parameterizing  the  sound  speed  as  a function  of  altitude: 


c m c ■ 300  m/sec  h <100  km 
o 


h >100  km 


(15) 
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so  that  the  sound  speed  Increases  linearly  to  2cq  at  an  altitude  pf 
160  km.  In  the  case  of  a sonic  boom  from  a Mach  2 aircraft,  the  upper 
turning  point  of  the  approximate  ray  trace  will  therefore  be  at  an 
altitude  of  160  km. 


The  pressure  and  density  in  our  model  atmosphere  may  again  be 
found  by  the  use  of  the  sound-speed  relation 


c2  « 1£ 

P 


and  the  hydrostatic  relation 


5h  P - -P8 


These  equations  immediately  lead  to  a local  scale  height 


H- 1 - - A S.  p - ^ 
s p h c2 


which  may  be  integrated  to  obtain  the  pressure  as  a function  of  altitude 
(Fig.  2b). 


E.  Sonic  Boom  Properties  as  a Function  of  Altitude 


Given  our  model  atmosphere  Eos.  (15)  and  (18)  and  our  propagation 
of  a sonic  boom  Eqs.  (3),  (A)  and  (7)  we  may  easily  numerically 
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integrate  to  find  the  sonic  boon  starting  at  an  altitude  of  17  Vcm,  with 

_3 

strength  sq  - 10  and  length  Lq  ■ 30  meters,  believing  these  to  be 
characteristic  of  a supersonic  aircraft.  (One  effect  not  mentioned 
previously  was  taken  into  account:  that  the  increase  of  sound  speed 
with  altitude  also  causes  a lengthening  of  the  boom,  contributing  to 
lowering  its  strength.) 

Figure  3 shows  the  strength  of  the  boom  as  a function  of 
altitude.  The  boom  reaches  a strength  of  little  more  than  0.1  at 
160  km  altitude.  Note  that  an  extrapolation  of  the  isothermal  atmosphere 
(valid  below  100  km)  would  have  .given  a strong  shock  in  this  region. 

Thus  the  thermospheric  properties  just  manage  to  avoid  the  production 
of  strong  shocks. 

From  this  result  we  may  calculate  thermospheric  heating  and  winds. 

Inside  the  shock  AT  — s T < 30° 

Y a 

2 

Behind  the  shock  AT  a ]L- ^ s3  T <0.03° 

a , , 3 a 

12y 

AT 

Wind  vS  ' 3 cm/sec 

a 

These  values  are  clearly  negligible  compared  to  other  fluctuations  in 
the  atmosphere. 
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Figure  4 shows  the  calculated  energy  in  the  sonic  boom  as  a 
function  of  altitude.  We  see  that  90%  of  the  energy  is  lost  before 
reaching  160  km  altitude. 

Figure  5 shows  the  lengthening  of  the  boom.  At  160  km  altitude 
the  boom  reaches  a length  of  about  4 km. 

We  may  describe  the  history  of  an  upward-going  sonic  boom,  following 
Figures  4 and  5,  by  approximating  the  true  refractive  path  (with  an 
upper  turning  point  at  about  160  km)  by  a vertically  upward  path  to 
160  km  followed  by  a vertically  downward  path  to  the  ground  from  160  km. 
In  that  approximation  only  3%  of  the  initial  energy  is  remaining  at  the 
upper  turning  point,  and  the  boom  has  lengthened  to  4 km.  One  can  show 
that  in  the  downward-going  leg  of  the  path,  only  about  30%  of  the 
energy  at  160  km  is  lost,  and  the  length  does  not  change  significantly; 
the  smallness  of  these  changes  is  due  to  the  rapid  decrease  in  the 
strength  as  the  boom  enters  denser  atmosphere.  Hence  the  signal  reach- 
ing the  ground  after  refracting  off  the  thermosphere  is  approximately 
50  times  less  energetic  than  the  sonic  boom  that  travelled  directly 
downward  from  the  aircraft,  and  it  is  about  4 km  in  length,  giving 
peak  frequencies  in  the  infrasonic  region  (-0.1  Hz). 

The  loss  in  energy  and  the  lengthening  results  in  a strength  at 
the  ground  of 
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If  this  signal  bounces  off  the  ground  and  undergoes  another  thermo- 
speric  refraction,  its  strength  will  remain  so  small  that  no  signif- 
icant energy  loss  or  lengthening  will  occur, 

F.  Other  Effects 


A more  extensive  treatment  of  this  problem  would  consider  the 
correct  direction-changing  path  of  an  angled  sonic  boom  that  spreads 
in  three  dimensions.  Such  a treatment  would  also  consider  the  caustic 
that  forms  after  the  wavefront  has  turned  over.  In  addition  viscosity 
in  the  air  could  be  included.  After  this  work  was  completed  an 
extensive  paper  by  Gardner  and  Rogers  (Naval  Research  Laboratory)  was 
received  that  took  all  these  effects  into  account.  They  reached 
conclusions  that  are  qualitatively  and  even  quantitatively  similar  to 
those  we  have  derived  with  our  simplified  model. 
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RESEARCH  NEEDS 


Widespread  interest  in  the  "East  Coast  Mystery  Booms"  has  focused 
attention  on  the  problem  of  understanding  the  propagation  of  low  frequency 
sound  in  a strongly  stratified,  inhomogeneous  atmosphere  having  winds  of 
various  strengths  at  different  levels.  In  the  United  States  the  subject 
received  considerable  attention  in  the  1950s  and  1960s,  first  because  of 
an  interest  in  detecting  nuclear  explosions  in  the  atmosphere  and  second 
because  of  concern  over  the  effects  of  flying  supersonic  transports,  both 
American  and  the  Concorde.  More  recently  the  subject  of  infrasound  has 
been  largely  neglected.  At  the  start  of  the  East  Coast  events  there  was 
only  one  low  frequency  acoustic  array  in  operation,  the  Lamont-Doherty 
facility.^  At  Lamont,  Donn  operates  two  tripartite  arrays  of  capacitor 
microphones.  The  larger  array  is  on  a 600  to  900  meter  spacing  and  has 
been  in  operation  since  1967  while  the  smaller  array  is  on  60  - 80  meter 
spacing  and  was  put  in  operation  in  July,  1977.  During  March  and  April, 
1978,  in  order  to  supplement  existing  observational  facilities,  Mitre  oper- 
ated three  two-microphone  correlators  at  Bedford,  Massachusetts,  Atlantic 
City,  New  Jersey  and  McLean,  Virginia.  The  study  of  the  East  Coast  Booms 
was  aided  by  the  existence  of  several  seismic  arrays  where  atmospheric 
pressure  disturbances  were  recorded  only  along  the  vertical  axis  of  the 
seismograph.  While  useful,  seismometers,  because  of  their  poor  perfor- 
mance in  the  0.01  to  5 Hz  region  and  because  of  the  Impedance  mismatch 


across  the  air-ground  interface,  are  less  than  ideal  Instruments  for 
observing  atmospheric  disturbances. 

The  recent  investigation  by  the  Lamont^’ and  Kiruna^ ’ ^ ^ groups 
have  indicated  that  numerous  artificial  sources  can  be  detected  at  large 
distances,  for  example  rocket  launches,  the  operations  of  SSTs,  paper 
mills,  hydroelectric  plants  and  offshore-oil  drilling  equipment.  The 
indicated  richness  of  these  sounds  suggest  that  a number  of  properties  of 
the  atmosphere,  for  example  stratospheric  and  thermospheric  winds,  could 
be  monitored  on  a routine  basis  using  netted  acoustic  arrays.  Such  obser- 
vations are  possible  because  of  a fundamental  property  of  the  atmosphere; 
infrasonic  waves  can  propagate  over  long  distances  because  the  attenuation 
is  very  small  in  the  0.01  to  5 Hz  region,  a frequency  band  where  natural 
atmospheric  noise  is  also  very  low. 

The  general  considerations  presented  above  suggest  a number  of  research 
priorities : 

A.  Observational 

1.  The  construction  and  operation  of  several  netted  infrasonic  arrays 

operating  in  the  0.01  to  5 Hz  region  to  monitor  Concorde  and  other 

artificial  sounds  of  infrasound  along  the  East  Coast.  The  existing 

Lamont  facility  needs  to  be  supplemented  by  other  arrays  in  order  to 

understand  better  the  Concorde  signals.  As  noted  by  Gardner  and 
9 

Rogers  an  array  on  a base  length  of  one  kilometer  has  poor  directional 
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resolution  of  sound  waves  in  the  few  tenths  of  a Hertz  range. 
Commercially  available  capacitor  microphones  are  relatively  inex- 
pensive, a few  hundred  dollars,  and  sophisticated  data  processing 
techniques  developed  in  other  fields  (radar,  ocean  acoustic)  could 
be  readily  applied  to  the  atmospheric  problem.  Both  fixed  and  por- 
table arrays  should  be  employed  to  determine  the  background  noise 
levels  in  the  infrasonic  frequency  region  as  a function  of  time  of 
day,  season  and  weather  as  well  as  array  location.  Such  data  are 
essential  if  we  are  to  understand  the  origin  of  the  natural  noise 
as  well  as  to  plan  the  location  of  future  arrays. 

2.  Source  spectra  of  artificial  sources.  The  near  field  spectra 
of  artificially  produced  infrasound  is  very  poorly  understood.  For 
example  the  spectra  of  conventional  jet  aircraft  below  60  Hz  is  not 
known  and  the  source  spectrum  in  the  infrasound  region  of  the  Concorde 
is  poorly  understood.  Data  on  the  infrasound  spectra  of  industrial 
facilities,  steam  power  plants  and  hydroelectric  facilities  are  not 
available.  These  spectra  are  needed  in  order  to  study  the  non-linear 
conversion  of  frequencies  in  a weak  shock,  the  attenuation  of  weak 
shocks  at  various  altitudes,  and  the  contribution  of  various  artificial 
sources  to  signals  in  the  infrasound  region  of  the  acoustic  spectrum. 

3.  Correlation  of  ionosphere  soundings  with  acoustic  detection. 

14 

Liszka  and  Olsson  report  that  flights  of  supersonic  fighters  (Saab 
35  Draken)  produced  both  infrasonic  (2  Hz)  signals  at  ground  level 
and  fluctuations  in  the  E layer  reflection  height  of  as  much  as  7 km. 
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The  interaction  of  weak  shocks  with  the  weakly  ionized  portion 
of  the  ionosphere  should  be  investigated  both  to  understand  the 
nature  of  the  coupling,  the  possible  contribution  of  the  neutral- 
ion  interaction  to  dissipation  at  high  altitudes  and  the  possible 
use  of  ionospheric  observation  to  detect  weak  shocks  launched  at 
ground  level  or  in  the  atmosphere.  Liszka  and  Olsson  also  observed 
a short-duration  b~ layer  echo  at  130  km  50  sec  before  the  arrival 
of  the  infrasonic  waves,  suggesting  that  the  supersonic  aircraft  may 
have  generated  another  mode  of  propagation,  possibly  gravity  waves. 

4.  Artificial  noise -generating  arrays.  An  array  of  infrasound 
sources  producing  a narrow  beam  of  high  intensity  could  provide  a 
powerful  tool  for  probing  the  physical  conditions  in  the  high  atmo- 
sphere. Such  a facility  would  be  particularly  useful  in  investigating 
the  dissipation  mechanisms  in  the  thermosphere.  Recent  calculations 

of  the  dissipation  assume  that  classical  viscosity  and  heat  conduction 

9 

are  the  principal  mechanisms  for  dissipation.  For  the  thermosphere, 
inhomegenieties  in  the  density  and  temperature  fields  on  temporal  and 
spatial  scales  smaller  than  those  of  the  transmitted  wave  can  act  as 
an  effective  eddy  viscosity  which  could  be  much  larger  than  the 
classical  viscosity.  A sound  array  coupled  to  an  ionosonde  facility 
would  be  a powerful  tool  for  investigating  the  neutral-ion  interac- 


tions noted  above 


Theoretical 


1.  Origin  of  background  noise  in  the  infrasound  region.  Natural 
infrasonic  noise  is  generally  attributed  to  turbulent  fluctuations 
in  the  wind  field.  The  relations  between  these  fluctuations  and 
ambient  weather  conditions  have  not  been  adequately  explored.  An 
understanding  of  the  nature  and  origin  of  the  natural  background 
noise  is  essential  to  selecting  quiet  sites  for  infrasonic  arrays. 

2.  Attenuation  of  weak  shocks  in  the  high  atmosphere.  The 
concepts  of  classical  viscosity  and  heat  conduction  are 
probably  no  longer  appropriate  when  the  collision  frequency  in  the 
ambient  atmosphere  is  of  the  same  order  as  the  infrasound  frequencies. 
The  physical  nature  of  the  attenuation  in  the  atmosphere  above  about 
100  km  needs  to  be  understood  in  order  to  interpret  long-range  pro- 
pagation of  infrasound. 

3.  Coupling  of  weak  shocks  with  ionized  particles.  The  observation 
of  E layer  height  variations  associated  with  the  upward  propagation 
of  a weak  shock  suggest  a strong  coupling.  The  detailed  physics  of 
the  coupling  needs  examination  as  does  the  possible  contribution  of 
weak  plasma  effects  to  the  attenuation. 

4.  Artificial  generation  of  gravity  waves  in  the  atmosphere.  The 
present  report  has  examined  the  infrasound  part  of  the  spectrum. 
Artificial  sources  of  gravity  waves  may  exist,  for  example  Concorde, 


and  these  artificial  waves  could  perturb  the  ionosphere.  A variety 
of  potential  artificial  sources  should  be  examined  in  terms  of  source 
strength  so  as  to  determine  whether  observations  below  0.01  Hz  could 
yield  information  about  artificial  sources  and  about  the  propagation 
of  gravity  waves  in  the  atmosphere. 
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SONIC  BOOM  PROPAGATION 


A.  Fundamental  Equations 

Consider  a sonic  boom  consisting  of  a shock  discontinuity  followed 
by  a continuous  waveform  of  length  L (Fig.  A-l).  The  shock  is  travelling 
into  ambient  atmosphere  at  speed  u: 


Pressure 


Figure  A-1  CONDITIONS  IN  A PROPAGATING  WEAK  SHOCK 

Just  after  the  shock  (i.e.  at  x small  but  positive)  we  define  the 

pressure,  density,  and  temperature  by  pQ  ■ p(o);  » p(o);  Tq  * T(o). 

The  velocity  of  the  gas  v(x)  is  non-zero  within  the  shock,  and  again  we 

define  v • v(o) . 
o 
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The  full  set  of  relations  between  quantities  on  either  side  of  the 


— — * 


shock  discontinuity  are  the  following: 


Mass  conservation 

pau 

' P (u  - V ) 
o o 

(1) 

Momentum  conservation 

p uv 

a 

“ P ~ P 
o ro  a 

(2) 

Energy  conservation 

P u 
a 

c (T  - T ) + 4 v2 
vo  a 2 o 

* p v 

0 0 

(3) 

Equation  of  state 

Po  - 

' PoVo 

(4) 

Within  the  sonic  boom  these  same  conservation  laws  yield: 


Pau  - P (x)  [ - v(x)]  (5) 

pau  axv0°  “ 5XP(X)  (6) 

p u 3 |c  (T  - T ) + 4 v2)  * 3 (pv)  (7) 

a x | v a l ] x 

p(x)  - p(x)  Rn  T(x)  (8) 


Equations  (1)  - (4)  yield  a number  of  relations  by  simple  algebra, 
which  we  present  in  terms  of  the  shock  strength 

si  (P0  - Pfl)/Pa  (9) 


> 


I 


O 


and  the  ambient  sound  speed 


~\ 


c - (YPa/Pa)1* 


(10) 
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One  can  show  chat  Che  shock  speed  is 


u.c(1+*is)‘ 


which  yields  for  s<<1 


The  temperature  just  behind  the  shock  is 


To-T,(lls)('  + 1#s)(1*1tH’1 


which  yields  for  s<<1 


AT/T  - ^ s 
a y 


The  entropy  change  f.cross  the  shock  is 


t 


AS  - cv  log  (To/Ta)  + ^ log  (oa/co) 


which  yields 


t 


t 


(ID 


(II-M 


(12) 


(11-1) 


(13) 
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B.  Relations  within  the  Sonic  Boom 

The  exact  manner  In  which  the  various  thermodynamic  variables  change 
as  a function  of  distance  behind  the  shock  is  a complicated  question  of 
fluid  dynamics.  However,  we  may  make  some  general  remarks: 

• The  integral  of  the  absolute  overpressure  must  ultimately 

be  close  to  zero  as  the  gas  is  assumed  to  return  to  velocity 
much  less  than  the  sound  speed  after  the  shock  has  passed. 

Thus  some  region  of  negative  overpressure  is  a necessity. 

• In  regions  of  smooth  pressure  behavior,  the  thermodynamic 
changes  will  be  isentropic.  However,  a second  shock 
behind  the  first  will  also  have  an  entropy  jump.  Because 
of  the  distortion  of  the  shape  of  the  shock,  the  entropy 
jumps  in  the  second  shock  will  be  less  than  that  of  the 
first  shock. 

We  must  determine  the  total  energy  in  a unit  area  of  wavefront,  and 
the  energy  left  behind  the  wavefront  in  the  ambient  atmosphere. 

The  energy  per  unit  volume  at  a distance  x behind  the  shock  is 


c = c (T  - T ) + \ v2  p(x)  (14) 

v a i 


* 


For  s«l  we  find 


C.  Relations  after  Sonic  Boom  Passage 


An  entropy  jump  across  the  first  shock  followed  by  lsentroplc  expan- 
sion back  to  ambient  pressure,  will  leave  the  atmosphere  heated  behind 
the  shock.  The  second  shock  in  an  ideal  N wave  will  generate  as  much 
entropy  as  the  first,  since  the  underpressure  in  the  following  shock  wave 
equals  the  overpressure  in  the  leading  wave.  In  an  actual  case,  the  shock 
lengthens  as  a result  of  momentum  transfer  and  the  second  shock  will  be 
weaker.  In  this  report  we  consider  the  weakening  of  the  shock  resulting 
only  from  entropy  jump  at  the  first  shock,  recognizing  this  approximation 
may  underestimate  the  fractional  energy  loss  by  as  much  as  a factor  of 
two  but  probably  much  less. 

AS  « cv  log  (T1/Ta)Y  (20) 


so  that 


V1. 


(1  + s),/v(ltIil.)(1+Iiis)-1 


(21) 


For  s<<l  we  have 


AT 

T 

a 


1 1=1 

12y  3 


0.03s3 


(I 1—2 ) 


•o  that  even  moderate-strength  shocks  heat  the  atmosphere  very  little. 
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The  energy  loss  In  travelling  a distance  L is  AE  * CyATp^L  , 
and  using  our  formulas  for  AT,  E , noting  that  cy  - ^/(y-1)  we  find 
the  fractional  energy  loss  in  traveling  a distance  L is 


AE  ~ _s 
E 2y 

Hence  we  may  immediately  write  the  energy  loss  equation 


(22) 


I - f - — 
E ?h  “*2yL 


(11-4) 


or 


Ii  psi- 
E h yL 


for  a N wave  maintaining  its  ideal  shape. 


The  momentum  transferred  to  the  atmosphere  is  related  to  the  energy 
deposition.  A wave  travelling  at  nearly  acoustic  speeds  (i.e.  s<<1)  would 
have  the  relation 

Energy  deposited 
Momentum  transferred  * c 

Thus  the  velocity  imported  to  the  atmosphere  after  the  passage  of  the 
boom  would  be 

c AT. 


Energy  deposited/unit  mass  m v a 
* c c 


(24) 
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This  lnmedlately  yields 
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